Diffusion of Nonequilibrium Quasiparticles in a Cuprate Superconductor 
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We report a transport study of nonequilibrium quasiparticles in a high-T c cuprate superconductor 
using the transient grating technique. Low-intensity laser excitation (at photon energy 1.5 eV) was 
used to introduce a spatially periodic density of quasiparticles into a high-quality untwinned single 
crystal of YBa2Cu306.5- Probing the evolution of the initial density through space and time yielded 
the quasiparticle diffusion coefficient, and both inelastic and elastic scattering rates. The technique 
reported here is potentially applicable to precision measurement of quasiparticle dynamics, not only 
in cuprate superconductors, but in other electronic systems as well. 

PACS numbers: 74.25.Gz, 78.47,+p 



Quasiparticles are the elementary excitations of a su- 
perconductor, created when a Cooper pair of electrons 
breaks apart. The dynamic properties of quasiparti- 
cles, that is their rates of diffusion, scattering, trap- 
ping, and recombination, are critical for applications 
of conventional superconductors in X-ray detectors QJ 
and in the manipulation of superconductor-based qubits 
0- In more exotic superconductors, such as the high- 
Tc cuprates, a better understanding of quasiparticle dy- 
namics may help to uncover the mechanism for Cooper 
pairing. A special property of the cuprate superconduc- 
tors is the d-wave symmetry of the gap function, which 
leads to an unusual quasiparticle spectrum. The mini- 
mum energy for the creation of a quasiparticle depends 
on the direction of its momentum 3J. It is zero for mo- 
menta in the 'nodal' direction, oriented at 45° relative to 
the Cu-0 bond. The most energetically expensive quasi- 
particles are the 'antinodaP ones, whose momenta are 
nearly parallel to the bond. The antinodal quasiparticles 
are the mystery particles of cuprate superconductivity. 
Because they feel the pairing interaction most strongly, 
their properties may hold the key to high-T c supercon- 
ductivity. Unfortunately, their tendency to form strong 
pairs makes them difficult to study. In thermal equilib- 
rium the population of quasiparticles is overwhelmingly 
dominated by the low energy nodal ones. As a result, 
transport measurements performed in equilibrium, such 
as microwave £j| and thermal [f| conductivity, are insen- 
sitive to antinodal quasiparticles. 

In this work a transient grating technique was de- 
veloped and used to probe the transport of nonequi- 
librium quasiparticles in the high-T c superconductor 
YBa2Cu3C>6.5. We find that the diffusion coefficient, 
D, is much smaller than the value obtained in measure- 
ments on equilibrium quasiparticles in the same material. 
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The disparity in D suggests that the current experiment 
probes a population of quasiparticles that are not near 
the nodes, and are perhaps close to the antinodal regions 
of momentum space. 

The nonequilibrium quasiparticles were introduced us- 
ing short optical pulses. To probe their propagation, we 
generated a spatially periodic population by interfering 
two pulses at the sample surface. The spatial period, X g , 
equals A/2sin(9, where A is the wavelength of the pulse 
and 29 is the angle between the two pump beams. The 
nonequilibrium quasiparticles cause a change in the index 
of refraction at the laser frequency |6j which is a linear 
function of their density pj. As a result, the sinusoidal 
variation in quasiparticle density creates an index grat- 
ing which can be detected by the diffraction of a probe 
pulse. 

After creation, the distribution of quasiparticles 
evolves due to the combined effects of recombination and 
diffusion. In the process of recombination a pair of quasi- 
particles jumps back into the Cooper pair condensate 
with the simultaneous transfer of their creation energy to 
some other form (e.g. phonons). The amplitude of the 
grating may also decay as quasiparticle diffusion drives 
the system towards a spatially homogeneous quasipar- 
ticle concentration. The goal of the experiment is to 
disentangle these effects and measure both the rates of 
recombination and diffusion. 

The transient grating technique has been used success- 
fully in a wide variety of applications, including exciton 
diffusion, dynamics of biomolecules, propagation of ultra- 
sound, and thermal diffusion 8]. Observing the propaga- 
tion of superconducting quasiparticles requires the ability 
to detect the transient grating at extremely dilute con- 
centrations. At the low excitation densities needed to 
detect their propagation, the quasiparticles produce no 
more than 1CP J fractional change in the index of refrac- 
tion. The diffracted intensity from such a grating would 
therefore be of order 10~ 10 of the incident probe intensity 
and consequently very difficult to detect. 



Detection of the grating required measurement of the 
amplitude of the diffracted wave, which is a part in 10 5 of 
the probe, rather than the intensity. This is accomplished 
through the use of heterodyne detection 0, [^} . The key 
element in the success of the heterodyne technique is the 
diffractive optic (DO) beamsplitter, which creates pairs 
of pump and probe beams [HI Il2| . The DO element in 
our experimental setup (Fig. 1) is an array of 10 sepa- 
rate 2.5 mm square phase masks, each with a different 
grating period, on a fused silica substrate. A beam from 
a TkSapphire laser (A=800 nm and pulse repetition rate 
80 MHz) is split into primary pump and probe beams, 
which are focused onto one of the phase masks at a small 
angle with respect to each other. The phase mask splits 
each of the primary beams by diffraction into the m=±l 
orders. (For clarity, only the probe beam paths are shown 
in Fig. 1). 

The interference of the two pump beams creates a spa- 
tially varying index of refraction in the sample which 
reaches a depth 1000 A below the surface and has half 
the period of the phase mask. We detect the index vari- 
ation using an implementation of the heterodyne tech- 
nique that enables absolute calibration of the phase 01 • 
Each of the two probe beams is specularly reflected from 
the surface of the sample, and returns, via the DO, to a 
Si detector. The two possible round-trip beam paths are 
shown as solid and dashed lines in Fig. 1. With the use of 
the DO beamsplitter, detection of the diffracted probe is 
automatically aligned. The diffracted component of PI is 
precisely colinear with reflected P2, and vice versa. The 
experiment is performed by alternately blocking one of 
the two reflected beams. If reflected PI is blocked, then 
reflected P2 and diffracted PI are mixed in the detector. 
This measurement by itself is insufficient to extract the 
wave amplitudes because the relative phase of PI and P2 
is undetermined. However, simply blocking reflected P2 
instead of PI produces a mixed signal with the conjugate 
phase. Comparing the detector output for the two con- 
jugate beam paths fixes the absolute phase and therefore 
the wave amplitudes as well. 

By adjusting the phase delay of PI relative to P2, the 
output of the Si photodiode measures either the change 
of the specular reflection coefficient due to the grating, 
R, or the amplitude of the diffraction efficiency, TG. R 
is proportional to the spatial average of the quasipar- 
ticle concentration, whereas TG is proportional to the 
component of the quasiparticle concentration at the fun- 
damental period of the grating, A g 0|. In Fig. 2 we plot 
both R and TG for X g =2 /im as a function of time delay, 
for several intensities of the excitation pulses, and conse- 
quently for a range of initial quasiparticle concentrations. 
The temperature of the sample was 5 K and the grating 
wavevector was oriented along the crystallographic b axis 
(parallel to the Cu-0 chains). 

The curves are normalized to their value at time de- 
lay zero. The R curves, which decay due to recombina- 
tion, are nonexponential, and their characteristic rate of 
decay increases with increasing concentration 0. The 
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FIG. 1: Illustration of the beam path for heterodyne tran- 
sient grating detection. Pump and probe beams from the 
laser are split at the diffractive optic (for clarity only the 
probe beams are shown). A spherical mirror and plane fold- 
ing mirror (represented schematically by a lens in the sketch) 
focus the beams to a single 100 fim spot on the sample. Af- 
ter specular reflection and diffraction at the sample surface, 
the two probe beams are recombined by the diffractive optic 
and directed to a Si photodiode detector. The wavevector of 
the quasiparticle density variation is changed, without opti- 
cal realignment, by translating the diffractive optic so that a 
different phase mask in the array is inserted in the beam. 



increase of recombination rate with density is consistent 
with the idea that each quasiparticle must encounter an- 
other to scatter into the Cooper pair condensate. The 
corresponding TG curves depend strongly on the pump 
intensity as well. However the time dependence at each 
pump intensity is different from R because TG reflects 
the combined effects of recombination and propagation. 
In systems where the recombination rate is independent 
of density, it is relatively straightforward to separate the 
effects of particle decay and diffusion. If the average con- 
centration decays exponentially with rate 7, the ampli- 
tude of the grating decays with rate 7 + Dq 2 , where D 
is the diffusion coefficient and q is the wavevector of the 
grating, 2ir/\ g . The ratio TG/R, which decays simply as 
exp(-Dq 2 t), isolates the effects of diffusion on the evolu- 
tion of the particle density. 

With this example in mind, we are led to consider 
TG/R for our data (see Fig. 3). The left, and right, 
panels show TG/R for grating periods of 2, and 5 /xm, 
respectively, for several values of the pump laser inten- 
sity. Unlike the example of density-independent recombi- 
nation, the decay of TG/R is nonexponential and highly 
dependent on the excitation density. At high intensity 
TG/R recovers its value at zero time delay, after an ini- 
tial rapid decrease. TG/R is nearly independent of grat- 
ing period at high intensity. As intensity decreases, the 
minimum of TG/R moves systematically towards longer 
times and the decay approaches a simple exponential. 
The exponential rate at low pump laser intensity clearly 
depends strongly on X g , as is expected if the grating de- 
cays due to diffusion. 

The dependence of TG/R on laser intensity and grat- 
ing period results from an interplay of density-dependent 
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FIG. 3: The ratio TG/R for the same pump intensities as 
in Fig. 2 (with P=0.63 mW omitted for clarity) for grating 
period 2/j, (left panel) and 5fi (right panel). At high excita- 
tion density the curves are nearly independent of the period. 
At low excitation density the grating enters the propagation- 
dominated regime where the decay of TG/R depends strongly 
on the period of the grating. 



FIG. 2: (top panel) Change in the specular reflection coeffi- 
cient, R, as a function of time delay following creation of the 
grating, for several values of the pump intensity. R is normal- 
ized to unity at time delay zero to illustrate the systematic 
slowing down of the recombination rate as the excitation den- 
sity is decreased, (bottom panel) Normalized amplitude of the 
diffracted probe beam as a function of time delay for the same 
values of pump intensity used to measure R. 



recombination, diffusion, and energy transfer. To unravel 
these effects we have modeled the quasiparticle dynamics 
by adding a quadratic recombination term to the diffu- 
sion equation: dn(x,i)/dt = Dd 2 n(x,t)/dx 2 — f3n 2 (x,t) 
|13|. Here, n is the quasiparticle density and (3 is the 
recombination coefficient, which is a measure of the in- 
elastic scattering rate [Tj|. The dynamics predicted by 
this equation depend on the relative magnitude of fin 
and Dq 2 . In the high density regime, (3n S> Dq 2 , re- 
combination dominates over diffusion. Because the rate 
is more rapid where n is larger, recombination distorts 
the sinusoidal grating by flattening the crests. The dis- 
tortion of the grating profile reduces the component of n 
at the fundamental grating period faster than its spatial 
average. This accounts for the initial decrease of TG/R. 
However, if this were the only process TG/R would never 
bounce back to its initial value, as it is seen to do in Fig. 
3. 

TG/R recovers in the high-intensity regime because 
the energy released as the quasiparticles recombine ap- 
pears in another form, one that also creates a change 
in the index of refraction. In conventional supercon- 
ductors the quasiparticle energy is converted to lattice 



vibrational energy, which (in the presence of electron- 
phonon coupling) will change the index. While phonons 
may play this role in the cuprate superconductors as well, 
excitations involving the flipping of electron spins may 
be generated instead when quasiparticles recombine [l6| . 
In either case, the recovery of TG/R follows if this sec- 
ondary form of energy does not diffuse on the 100 ps 
time scale. In this case, the original sinusoidal distribu- 
tion of energy is written into a stationary form of energy 
before the quasiparticles have had a chance to diffuse. 
When almost all the energy has been handed off from 
the quasiparticles to non-propagating modes, the grating 
recovers its sinusoidal profile and TG/R returns to its 
value at time zero. The energy remains frozen until the 
nanosecond time scale, when thermal diffusion becomes 
significant |l7j . 

The effects of diffusion dominate in the low density 
regime, f3n <C Dq 2 , where quasiparticle motion washes 
out the grating before energy transfer can take place. Ac- 
cording to the equation governing n(x, t), TG/R becomes 
q dependent and intensity independent in this limit. To 
find the rate of diffusion, we measured TG/R at low 
power (P=0.32 mW) for several grating periods between 
2 and 5 /im. The initial decay rate of TG/R at low power 
is plotted as a function of q 2 in Fig. 4. The decay rate de- 
pends systematically on the grating period, demonstrat- 
ing that the dynamics are in the propagation-dominated 
regime. The rates for the grating oriented along both the 
a and b crystallographic directions are plotted. For both 
directions the rate is a linear function of q 2 , demonstrat- 
ing that the quasiparticle propagation is diffusive. From 
the slope of a linear fit we determine that _D a =20 cm 2 /s 



4 




q 2 (10 8 cm" 2 ) 



FIG. 4: Initial decay rate of TG/R, for the same low intensity 
as in Fig. 3, plotted as a function of the square of the grating 
wavevector. Results are shown for two perpendicular orienta- 
tions of the grating. The linear dependence of the rate on q 2 
demonstrates that the propagation is diffusive. The slopes of 
linear fits yield diffusion coefficients D a = 20 cm 2 /s and Dt, 
= 24 cm 2 /s. 



and D},=2A cm 2 /s. The intercept as q tends to zero is 
the decay rate due to recombination. 

It is possible to infer the mean free time r and mean 
square velocity < v 2 > of nonequilibrium quasiparticles 
from the values of D quoted above. A lower bound on r 
is obtained by inserting the maximum quasiparticle ve- 
locity, the Fermi velocity vp, into the kinetic formula, 
D =< v 2 > t/2. The literature value p| Vp = 2 x 10 7 
cm/s yields a lower bound of 100 fs. An upper bound 
on t of essentially the same value can also be inferred 
from the experiment. If r were significantly longer than 
100 fs then quasiparticle propagation would be ballistic 
on the subpicosecond time scale. However, the time and 
wavevector dependence of TG proves that quasiparticle 
motion is diffusive at the earliest times we can resolve, 
which is ~ 300 fs after creation of the grating. Thus the 
allowed values of r and v are narrowly bracketed near 100 
fs and Vp, respectively. Furthermore, the measurements 
indicate that r is determined by elastic, rather than in- 



elastic, processes. The nonequilibrium particles survive 
for « 100 ps, and therefore scatter « 1000 times before 
decaying. This would be impossible if each scattering 
event resulted in a significant reduction of the quasipar- 
ticle's energy. 

Equilibrium measurements find (in the same crystal 
and at the same temperature) that r = 20 ps [l9j , which 
is 200 times longer than r of the nonequilibrium quasi- 
particles. The contrast suggests that the nonequilibrium 
quasiparticles are different from those present in ther- 
mal equilibrium. The equilibrium particles occupy states 
which have a small energy (fc^T) relative to the chemi- 
cal potential, and therefore lie very near the gap nodes. 
If the nonequilibrium quasiparticles are different, they 
must occupy higher energy states, perhaps ones closer to 
the antinodal regions of momentum space. It is possible 
that the constraints of momentum and energy conserva- 
tion prevent relaxation of antinodal quasiparticles to the 
nodal regions. If the particles are antinodal, it is relevant 
to compare the width of the antinodal quasiparticle peak 
as measured by photoemission with K/t as measured by 
nonequilibrium transport. The peak width of 14 meV 
[2(i| corresponds to t=50 fs, which is close to the trans- 
port value, particularly considering that photoemission is 
performed on Bi2Sr2CaCu20s rather than YBa2Cu306.5- 
The similarity of lifetimes suggests that the peak width 
in photoemission may be controlled by elastic scattering 
as well. 

The transient grating method reported here promises 
to be broadly applicable to superconductors, as well as 
other materials in which there is a gap in the quasiparticle 
spectrum. The technique works readily in transmission 
or reflection geometry and therefore can be applied to 
bulk materials or thin films. The propagation of quasi- 
particles can be tracked in any system where nonequi- 
librium excitations generate a change in the index of 
refraction at the laser wavelength. In conventional su- 
perconductors, quasiparticle diffusion can be measured 
without fabricating trapping layers and junction detec- 
tors. In more exotic systems with multiple or anistropic 
gaps, such as reported in this paper, the transient grat- 
ing technique can track the propagation of quasiparticles 
that conventional transport methods cannot detect. 

This work was supported by DOE-DE-AC03- 
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